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Kaon Condensation in a Nambu— Jona-Lasinio (NJL) Model at High Density 
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We demonstrate a fully self-consistent microscopic realization of a kaon-condensed colour-flavour 
locked state (CFLK") within the context of a mean-field NJL model at high density. The properties 
of this state are shown to be consistent with the QCD low-energy effective theory once the proper 
gauge neutrality conditions are satisfied, and a simple matching procedure is used to compute the 
pion decay constant, which agrees with the perturbative QCD result. The NJL model is used to 
compare the energies of the CFLK" state to the parity even CFL state, and to determine locations 
of the metal/insulator transition to a phase with gapless fermionic excitations in the presence of a 
non-zero hypercharge chemical potential and a non-zero strange quark mass. The transition points 
are compared with results derived previously via effective theories and with partially self-consistent 
NJL calculations. We find that the qualitative physics does not change, but that the transitions are 
slightly lower. 



I. INTRODUCTION 

Recently there has been interest in the structure of 
matter at extremely high densities, such as might be 
found in the cores of neutron stars. At large enough den- 
sities, the nucleons are crushed together and the quarks 
become the relevant degrees of freedom. The asymptotic 
freedom of QCD ensures that the theory is weakly cou- 
pled at high enough densities. This allows one to per- 
form weak-coupling calculations at asymptotically high 
densities. Such calculations have established that the 
structure of the ground state of quark matter is a colour 
superconductor (see for example Q, S S H 1 i, H i, S 
EPj llj J^i J4a J^J)- In particular, at densities high 
enough that the three lightest quarks can be treated as 
massless, the ground state is the colour-flavour-locked 
(CFL) state in which all three colours and all three 
flavours participate in maximally (anti)-symmetric pair- 

ing i3|i3,lil|i3. 

Determination of the QCD phase structure at moder- 
ate densities and in the presence of non-zero quark masses 
has proceeded in several ways. One approach has been to 
formulate a chain of effective theories, and then to match 
coefficients across several energy scales through these ef- 
fective theories to perturbative calculations. Coefhcients 
in the low-energy chiral effective theory are matched 
to calculationsperformed in high-density effective theo- 
ries (HDET) li [1II23 which in turn are matched to 
weakly-coupled QCD. This allows one to determine the 
properties of the Goldstone bosons and determine the ef- 
fects of small quark masses jMIHIllllllllllillllll 
|2^|33- Within this framework, it has been noted that, 
in the presence of a finite strange quark mass, neutral 
"kaons" (the lightest pseudo- Goldstone modes at high 
density with the same quantum numbers as their vacuum 
counterpart) can Bose-condense in the CFL state to form 
a kaon-condensed CFLK" phase with lower condensation 
energy PHHll. 

Unfortunately, the low-energy effective theory is only 
reliable for small perturbations and at moderate densi- 
ties the strange quark mass is not a small perturbation. 



A recent attempt has been made to extrapolate to large 
strange quark mass (mg) [s^l, but this approach has not 
dealt with additional complications in the condensate 
structure that allow different gap parameters for each 
pair of quarks. 

To deal with moderate quark masses, another ap- 
proach has been to study Nambu-Jona-Lasinio (NJL) 
models |3^ of free quarks with contact interactions 
that model instanton interactions or single gluon ex- 
change. These models are amenable to a mean-field treat- 
ment and exhibit a similar symmetry breakingpattern to 
QCD which results in CFL ground states 0,0. 

Within these models, one can study the effects of 
moderate quark masses through self-consistent solutions 
of the mean-field gap equations. This has led to a 
plethora of phases. In particular, several analyses show 
a transition to a colour-flavour locked phase with gapless 
fermionic excitations (the gCFL phase). These include 
both NJL-based calculatio ns .3L. 38b. .39, .40.) and effective- 
theory-based calculations |34 l4l| . Until recently, how- 
ever, the NJL calculations have excluded the possibil- 
ity of kaon condensation (see however 42] which consid- 
ers kaon condensation in the NJL model at low density), 
while the effective theories do not consider the compli- 
cated patterns in which the condensate parameters evolve 
at flnite quark masses. 

The goal of this paper is to show that one can com- 
bine the analysis of the low-energy effective theories, 
which exhibit kaon condensation, with the self-consistent 
mean-fleld analysis of the NJL model, which accounts 
for the full condensate structure. In particular, we use 
an NJL model based on single gluon exchange to find 
self-consistent solutions that correspond to the CFLK" 
phase; we show that these phases agree with the predic- 
tions of the low-energy effective theory; and we deter- 
mine how and where the zero temperature phase transi- 
tion to a gapless CFL phase occurs as one increases the 
strange quark mass. In addition, unlike previous work 
on the NJL model, our numerical solutions are fully self- 
consistent: we include all condensates and self-energy 
corrections required to close the gap equations. 
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We first describe the pattern of symmetry breaking 
that leads to the CFL and CFLK° states (Section HJl . 
Then we present our numerical results, demonstrating 
some properties of these states and determining the lo- 
cations of the zero-temperature phase transitions (Sec- 
tion After a careful description of our model (Sec- 
tion lIV|l we derive the low-energy effective theory, paying 
particular attention to the differences between QCD and 
the NJL model (Section 0). Here we demonstrate that, 
for small perturbations, our numerical solutions are well 
described by the effective theory, and we use our numer- 
ical results to compute the pion decay constant which 
agrees with the perturbative QCD results. Specific nu- 
merical details about our calculations and a full descrip- 
tion of our self-consistent parametrization are given in 
Appendix fXl 

We leave for future work the consideration of finite 
temperature effects, the analysis of the gapless CFLK" 
(gCFLK*^), the inclusion of instanton effects, the inclu- 
sion of up and down quark mass effects, and the possi- 
bility of other forms of meson condensation. 



II. COLOUR FLAVOUR LOCKING (CFL) 

QCD has a continuous symmetry group of U(1)b (Xi 
SU(3)l O SU(3)r (g) SU(3)c. In addition, there is an ap- 
proximate U(1)a axial flavour symmetry that is explicitly 
broken by anomalies. At sufficiently high densities, how- 
ever, the instanton density is suppressed and this sym- 
metry is approximately restored. 

The CFL ground state spontaneously breaks these con- 
tinuous symmetries through the formation of a diquark 
condensate 

{^Ll5^^pb) « ^se'^^'eabk + AeiSy^, + 5^5^). (1) 

The symmetry breaking pattern (including the restored 
axial U(1)a symmetry) is thus^ 

U(3)l0U(3)r0SU(3)c ^„ ^„ 

> bU(3)L+R+C®^2<Xl^2 (^) 

^3 

where the Z2 symmetries correspond to ipL —i^L and 
—ipR- It has been noted that the symmetry break- 
ing pattern at high density ij^Jl is the same as that that for 
hyper- nuclear matter at low density 43] . This leads one 
to identify the low-energy pseudo-scalar degrees of free- 
dom in both theories. We shall refer to the pseudo-scalar 
Goldstone bosons in the high-density phase as "pions" 
and "kaons" etc. when they have the same flavour quan- 
tum numbers as the corresponding low-density particles. 



The Z3 factor mods out the common centres. See ilOt for the 
exphcit representation. 



The CFL state preserves parity, and is preferred 
when instanton effects are considered. Excluding instan- 
ton effects, there is an uncountable degeneracy of phys- 
ically equivalent CFL ground states that violate parity. 
These are generated from the parity even CFL by the 
broken symmetry generators. 

The symmetry breaking pattern Q breaks 18 genera- 
tors. The quarks, however, are coupled to the eight glu- 
ons associated with the SU(3)c colour symmetry and to 
the photon of the U(1)em electromagnetism (which is a 
subgroup of the vector flavour symmetry) . Eight of these 
gauge bosons acquire a mass through the Higgs mech- 
anism and the coloured excitations are lifted from the 
low-energy spectrum. There remain 10 massless Nambu- 
Goldstone excitations: a pseudo-scalar axial flavour octet 
of mesons, a scalar superfluid boson associated with the 
broken U(1)b baryon number generator, and a pseudo- 
scalar 7y' boson associated with broken axial U(1)a gen- 
erator. There remains one massless gauge boson that 
is a mixture of the original photon and one of the glu- 
ons 0, 01 ■ With respect to this "rotated electromag- 
netism" U(1)q the CFL state remains neutral [45j| . 

The degeneracy of the vacuum manifold is lifted by the 
inclusion of a non-zero strange quark mass mj. In the 
absence of instanton effects and other quark masses, the 
ground state is not near to the parity even CFL state iQ), 
but rather, is a kaon rotated state CFLK". As — > 
this state approaches a state on the vacuum manifold 
that is a pure kaon rotation of the parity even CFL 

Even in the absence of quark masses, the vacuum 
manifold degeneracy is partially lifted by the anomalous 
breaking of the U(1)a axial symmetry which we have ne- 
glected: Instanton effects tend to disfavour kaon conden- 
sation by favouring parity even states, and thus delay the 
onset of the CFLK° until reaches a critical value (pos- 
sibly excluding it). The effects of anomaly and instanton 
contributions have been well studied 
and play an important quantitative role in the phase 
structure of QCD. Non-zero up and down quark masses 
also tend to disfavour kaon condensation. 

For the purposes of this paper, we shall neglect both 
the effects of instantons, and the effects of finite up and 
down quark masses. This will ensure that kaon condensa- 
tion occurs for arbitrarily small uis- Both of these effects 
open the possibility of a much richer phase structure, 
including condensation of other mesons (see for exam- 
ple j32i,.4£j). Future analyses should take these numeri- 
cally important effects into account, both in the effective 
theory and in the NJL model. 

The primary source of for kaon condensation is the fi- 
nite strange quark mass. To lowest-order, this behaves as 
a chemical potential l^ll^,!!^ (see © and (0J). In this 
paper, we also consider the addition of a hypercharge 
chemical potential as this removes many complications 
associated with masses and leads to a very clean demon- 
stration of kaon condensation. 
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III. SELF-CONSISTENT SOLUTIONS 

We consider four qualitatively different phases: Two 
are self-consistent mean-field solutions to the NJL model 
with a finite hypercharge chemical potential parameter 
/iy; the other two are self-consistent mean- field solutions 
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to the NJL model with a finite strange-quark mass pa- 
rameter TTis. In each of these cases, one solution corre- 
sponds to a parity even CFL phase and the other cor- 
responds to a kaon-condensed CFLK° phase. Our nor- 
malizations and a complete description of the model are 
presented in Section IIVI A full description of all the 
parameters required to describe these phases along with 
some typical values is presented in Appendix IXI 
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TABLE L Leading order shifts in the chemical potentials 
of the various quarks in the CFL and CFLK° states in the 
presence of a hypercharge chemical potential shift py. This 
follows directly from 15UII . 



I • 1 

450 500 550 

\p\ (MeV) fiY =0 
80 I ' 1 




450 500 550 

IpI (MeV) ^IY = My 

FIG. 1: Lowest lying quasiparticle dispersions about the 
Fermi momentum pp ~ i-iq = 500 MeV for the CFL phase 
with different values of the hypercharge chemical. All disper- 
sions have left-right degeneracy: we now consider the colour- 
flavour degeneracy. In the top plot fiy ~ 0, and the lowest 
dispersion has an eight-fold degeneracy and a gap of Ao — 25 
MeV. The upper band contains a single quasiparticle pairing 
(ru,gd,bs) with a gap of 4A6 -I- 2A3 = 54 MeV. In the middle 
plot, fiY = l^vf^ = 12.5 MeV, and (rs,bu) and (gs,bd) pairs 
are shifting as indicated in Table |I] In the last plot, two pairs 
have become gapless marking the CFL/gCFL transition. 
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IpI (MeV) fiY = 1.20My 

FIG. 2: Lowest lying quasiparticle dispersions about the 
Fermi momentum pp = ^J.q = 500 MeV for the CFLK" 
phase with different values of the hypercharge chemical. (The 
/iy = dispersions are the same as in the top of Figure 0) 
Again, all dispersion have a left-right degeneracy. In the top 
plot at /iy = ^iy /2 = 12.5 MeV, the eight-fold degenerate low- 
est band has split into eight independent dispersions. To lead- 
ing order in the perturbation, the splitting is described by Ta- 
ble Q but the lack of degeneracy indicates that there are also 
higher order effects. The lower plot at at ~ 1.20/^y « 30 
MeV is close to the CFLK'^/gCFLK'^ transition. The gapless 
band now contains only a single mode and is charged. 
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FIG. 3: Physical gap of the lowest lying excitation as a func- 
tion of the hypercharge chemical potential. The dotted line 
corresponds to the CFL phase: the phase transition to the 
gCFL occurs at fiy ~ I-i-y where the gap vanishes. The solid 
line corresponds to the CFLK" state. The transition to a 
gapless phase is delayed by a factor of 1.22. 

A. Finite Hypercharge Chemical Potential 

The CFL phase in the presence of a hypercharge chem- 
ical potential corresponds to the fully gapped CFL phase 
discussed in (s^- Here one models the effects of the 
strange quark through its shift on the Fermi surface 
Pf ~ ^J'q of the strange quarks. This can be seen by 
expanding the free-quark dispersion 



V^+Ml^\p\ + -^ + --. (3) 

or, more carefully, by integrating out the antiparticles 
to formulate the High-Density Effective Theory. (See for 
example 9, 19, 20, 30].) These leading order effects are 
equivalent to adding a hypercharge chemical potential of 
magnitude 

f^Y = (4) 
and a baryon chemical potential shift of 

6^iB = -. (5) 

We consider only the effect of the hypercharge chemi- 
cal potential here, holding ng fixed. Note that the rel- 
evant parameters are and iiq rather than and 
fJ-s = (J'b/'^- Ms is the constituent quark mass that 
appears in the dispersion relation whereas rn^ is the 
bare quark mass parameter; likewise, fiq is the corrected 
quark chemical potential that determines the Fermi sur- 
face whereas /is = /J_b/3 is the bare baryon chemical po- 
tential. (These distinctions are important because our 
model takes into account self-energy corrections.) 

The CFL phase responds in a trivial manner to a 
hypercharge chemical potential: the quasiparticle dis- 



FIG. 4: Physical gap of the lowest lying excitation as a func- 
tion of the strange quark mass. The dotted line corresponds to 
the CFL phase and the solid line corresponds to the CFLK" 
phase. We have normalized the axes in terms of = Ao 
for comparison with the hypercharge chemical potential case. 
The CFL/gCFL transition occurs at a slightly smaller value of 
M^Vm? « 45.4 MeV than the value of 46.8 MeV in IsTllalliot. 
This is due to the effects of the other parameters on the quasi- 
particle dispersion relations. We note that, as with fiy, the 
transition from the CFLK'' to a gapless phase is delayed rel- 
ative to the CFL/gCFL transition, but by a slightly reduced 
factor of 1.2. This is in qualitative agreement b ut q uantita- 
tive disagreement with the factor of 4/3 found in |34| . The is 
most-likely the result of our fully self-consistent treatment of 
the condensate parameters. 

persions shift such that the physical gap in the spec- 
trum becomes smaller; none of the other physical prop- 
erties change. In particular, as the hypercharge chem- 
ical potential increases, the coloured chemical potential 
/ig — — /iy decreases to maintain neutrality. The values 
of all of the gap parameters, the self-energy corrections, 
the densities and the thermodynamic potential remain 
unchanged until the physical gap in the spectrum van- 
ishes. (The apparent change in the magnitude of the gap 
parameters in the first figure of 0| is due to the shift 
in the baryon chemical lO which occurs if one uses the 
strange quark chemical potential shift rather than a 
hypercharge shift /iy.) This is a consequence of the Q 
neutrality of the CFL state js^]. In particular, the elec- 
tric chemical potential remains zero fie ~ and the state 
remains an insulator until the onset of the gapless modes. 
The same phenomena has also been noticed in the two- 
flavour case 0, Ell' As such, we can analytically 
identify the phase transition to the gCFL phase which 
occurs for the critical chemical potential 

My = Ao (6) 

where Aq — A3 — Ag is the physical gap in the spectrum 
in the absence of any perturbations. Throughout this 
paper we use parameters arbitrarily chosen so that /iy — 
Aq = 25 MeV to correspond with the parameter values 
in [23, HE 1^ • We show typical quasiparticle dispersion 
relations for this state in Figure ^ 
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IpI (MeV) M|/(2^,) = 0.83ai?, 

FIG. 5: Lowest lying quasiparticle dispersion relationships 
about the Fermi momentum pp ~ ~ 500 MeV for the CFL 
phase with two different values of the strange quark mass. 
(The Ms = dispersions are the same as in the top of Fig- 
ure Qualitatively this has the same structure as Figure Q 
except that middle dispersion is now split by higher order 
mass effects. 



The splitting of the dispersions can also be easily un- 
derstood from the charge neutrality condition ifSUI) and 
the leading order effects are summarized in Table Af- 
ter setting fig = —fJ.Y, the chemical potentials for the rs 
and gs quarks shift by —fiy whereas for the bu and bd 
quarks it shifts by +hy- Thus, the (gs,bd) and (rs,bu) 
pairs are the first to become gapless. 

The kaon-condensed hypercharge state is more com- 
plicated. One can again use the appropriate charge neu- 
trality conditions H50|l to estimate how the quarks will be 
affected by fiy, but the naive results hold only to lowest- 
order. In particular, the condensates of the CFLK" state 
also vary as /iy increases (see Table 1111}) . These higher 
order effects break all the degeneracy between the quark 
species and Figure [21 has nine independent dispersions. 

We shall compare the thermodynamic potentials of 
these two states later (see Figures [7| and |Hl, but we 
point out here that the transition to a gapless colour- 
flavour-locked state with kaon condensation (gCFLK°) 
occurs at a larger hypercharge chemical potential than 
the CFL/gCFL transition. This can be most easily seen 
in Figure 13 This is in qualitative agreement with 
and [IJ], but in quantitative disagreement. 




IpI (MeV) M|/(2/i,) = 0.84 ii'y 

FIG. 6: Lowest lying quasiparticle dispersion relationships 
about the Fermi momentum pp ~ i-iq = 500 MeV for the 
CFLK" phase with two different values of the strange quark 
mass. (The Ms = dispersions are the same as in the top 
of Figure Qualitatively this has the same structure as 
Figure H 



In the CFL/gCFL transition, two modes become gap- 
less simultaneously: the lower branches of the (rs,bu) and 
(gs,bd) pairs. One of these modes is electrically neutral 
(gSjbd) and it crosses the zero-energy axis giving rise to 
a "breach" in the spectrum. The other mode is electri- 
cally charged: as soon as in crosses, the electric chemical 
potential must rise to enforce neutrality. The state now 
contains gapless charged excitations and becomes a con- 
ductor. The result is that the the neutral gapless mode 
has two linear dispersions while the charged gapless mode 
has a virtually quadratic dispersion when electric neu- 
trality is enforced. (This was discovered in |37| and is 
explained in detail in 

In the CFLK°/gCFLK° transition, a single charged 
mode becomes gapless.^ Thus, immediately beyond the 



^ This mode pairs rs, gu, and bu quarks in quite a non-trivial man- 
ner. In the CFL, the quasiparticles form a nice block-diagonal 
structure in which the quarks exhibit definite pairing between 
two species. In the CFLKO, the block structure is more com- 
plicated and the pairing cannot be simply described: the lowest 
lying quasiparticle is a linear combination of the three rs, gu, 
and bu quark. 
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transition, the corresponding gCFLK" state will also be a 
conductor but there will be a single charged gapless mode 
with almost quadratic dispersion. Additional modes will 
continue to lower until either more modes become gap- 
less, or a first order phase transition to a competing phase 
occurs. 



B. Finite Strange Quark Mass 

The second pair of CFL/CFLK° states that we con- 
sider are self-consistent solutions to the gap equation in 
the presence of a finite strange quark mass. Qualitatively 
we expect to see similar features to the states at finite hy- 
percharge chemical potential and indeed we do as shown 
in Figures |31 and El 

Quantitatively, we notice a few differences with previ- 
ous analyses concerning the locations of the phase transi- 
tions to gapless states. Our parameters have been chosen 
to match the parameters in [s^, IH, ^3 ■ They find that 
the gCFL/CFL transition occurs at = 46.8 MeV, 

but the CFL/gCFL transition happens noticeably ear- 
lier with our model at M^/m = 43.9 MeV. This is due 
to a corresponding six-percent reduction in the conden- 
sate parameters and represents the effects of performing 
a fully self-consistent calculation. 

Another difference concerns the appearance of gapless 
modes in the CFLK° state (see Fig. 0J. This transition 
occurs at M^/fi = 52.5 MeV in our model — a factor of 
1.2 larger than the CFL/gCFL transition. This is some 
ten percent smaller than the factor of 4/3 derived in [s^ . 
This is likely due to the more complicated condensate 
structure we consider and the inclusion of self-energy cor- 
rections. 



IV. NJL MODEL 

We base our analysis on the following Hamiltonian for 
the NJL model 



J^^l (a • p - /X + 7oM) ^pfi + Hi, 



(7) 



Here we consider 9 species of quarks = 3 colours x 3 
flavours: Including the relativistic structure, there are 36 
quark operators in the vector ip. The matrices fi and 
M are the quark chemical potentials and masses respec- 
tively. 

We take the interaction to be a four-fermion contact 
interaction with the quantum numbers of single gluon 
exchange:'^ 



^^int — . 



(V'7^A'^V)('/'7 



(8) 



^ Here the matrices are the eight 3x3 Gell-Mann matrices and 
the -)''' are the Dirac matrices which we take in the chiral basis. 



The Gell-Mann matrices act on the colour space and the 
flavour structure is diagonal. We point out that this form 
of NJL interaction has the desirable feature of explic- 
itly breaking the independent colour SU(3)cl left and 
SU(3)cR right symmetries that some NJL models pre- 
serve. This is important because the condensation pat- 
tern does not explicitly link left and right particles: 
Our model thus has the same continuous symmetries as 
QCD, and the only complication to deal with is the gaug- 
ing of the single colour SU(3)c symmetry. 

Our goal here is to provide a non-perturbative model to 
discuss the qualitative features of QCD at finite densities. 
We model the finite density by working in the grand ther- 
modynamic ensemble and introducing a baryon chemical 
potential for all of the quarks: 



MS, 



(9) 



With only this chemical potential and no quark masses, 
our model has an U(3)L®U(3)R(g)SU(3)c/Z3 continuous 
global symmetry in which the left-handed quarks trans- 
form as (3, 1, 3) and the right handed quarks transform 
as (1, 3, 3). In the chiral basis we have explicitly'* 









(10) 



where F and C are SU(3) matrices. For an attractive 
interaction, this NJL model exhibits the same symme- 
try breaking pattern as QCD Q with a restored ax- 
ial symmetry. The difference between this NJL model 
and QCD is that the NJL model contains no gauge 
bosons. Thus, there are 18 broken generators which cor- 
respond to massless Goldstone bosons, and none of these 
is eaten. To effectively model QCD, we must remove 
the extra coloured Goldstone bosons. At the mean- field 
level, this is done by imposing gauge neutrality condi- 
tions 13 [S^. Once the appropriate chemical po- 
tentials are introduced, the dependence on the vacuum 
expectation values of the coloured Goldstone modes is 
cancelled and the low-energy physics of the NJL model 
matches that of QCD. 

The usual NJL model has a local interaction, but this 
is not renormalizable and needs regulation. For the pur- 
poses of this paper, we introduce a hard cutoff on each 



Our normalizations and conventions are 



TrfA^A 



1 cAB 
2^ ' 



75 = *70')'lT2T3 



1 



We also use natural units where c = h = ks = 1- 
* From this explicit representation we can see how the centres of 
the colour and flavours overlap giving rise to the Z3 factor. 
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of the momenta = 6'(A — |li5||) to mimic the effects of 
asymptotic freedom at large momenta: 



(2^)- 



J d^pd^p'd^qd^q' ApAp,A,Ag,x 



X <5(3)(p - p' + q - q')(Vii57''A^Vp')(Viq7^A^V^')- 

To study this model we perform a variational calculation 
by introducing the quadratic Hamiltonian 



Ho 



where 



d^p 



(27r)3 



3V4 



^lE(p)^j5 + BV;j5 + V-^bV; j (11 



E(p) = (3 • p - /X + - A 



(12) 



and then computing the following upper bound 56] on 
the thermodynamic potential Q of the full system: 



(13) 



f2o is the thermodynamic potential of the quadratic 
model and the expectation value ()q is the thermal av- 
erage with respect to the quadratic ensemble defined by 
Hq. In principle, the quadratic model is exactly solvable, 
thus the upper bound can be computed. One then varies 
the parameters A and B to minimize this upper bound, 
obtaining a variational approximation for the true en- 
semble. In the zero-temperature limit, this is equivalent 
to simply minimizing the expectation value of the Hamil- 
tonian over the set of all Gaussian states. 

In practice, it is difficult to vary with respect to all 
possible quadratic models since the space is of uncount- 
able dimensionality. In this paper we restrict ourselves 
to minimizing over homogeneous and isotropic systems. 
This is equivalent to performing a fully self-consistent 
mean- field analysis. The condition for the right hand 
side of (|13|l to be stationary with respect to the varia- 
tional parameters is equivalent to the self-consistent gap 
equation. 

The microscopic analysis presented in this paper con- 
sists of choosing reasonable parameterizations of A 
(which includes the chemical potentials, masses and re- 
lated corrections) and B (which includes the gap param- 
eters A) that are closed under the self-consistency con- 
dition, and numerically finding stationary points of this 
system of equations. (As A and B are arbitrary 36 x 36 
matrices subject only to A = A^ and B = — B^, a full 
parametrization consists of 2556 parameters and was too 
costly for the present author to consider. However, the 
parametrization chosen is quite natural and fully closed.) 
Once the parameters A and B are found, the properties 
of the ensemble can be computed by diagonalizing the 
quadratic Hamiltonian. 

As discussed in Section IVEl and [isf. we must impose 
the appropriate gauge charge neutrality conditions. This 
is done by introducing bare gauge chemical potentials 
into the model and choosing these to ensure the final 
solution is neutral. 



To impose a charge neutrality condition, we instead 
vary f^n (along with with the other parameters) to obtain 
a neutral solution (again we note that the total charge 
and other correlations of the state depend only on the 
corrected parameters fin). Once this solution is found, Sfi 
is computed and the required bare chemical potential fi = 
fiR — Sfj, determined. Despite the fact that the self-energy 
corrections depend only on the corrected parameters {fin 
etc.), the thermodynamic potential depends on both the 
corrected and the bare parameters and so this last step 
is important. 

One must also be careful about which thermodynamic 
potential is used to compare states when neutrality con- 
ditions are enforced as we are no longer in the grand 
ensemble. The differences between the potentials of the 
relevant ensembles are proportional to terms of the form 
Qfi, however, so for neutrality conditions, Q = 0, and the 
thermodynamic potential may still be used to compare 
states. 



A. Numerical Techniques 

We sketch here the method used to calculate the ther- 
modynamic potential and perform the variational min- 
imization. First, we express the Hamiltonian H in the 
following simplified form 



(14) 



where H is a hermitian matrix. In order to do this and 
include the "anomalous" correlations (ip'^P) , we must use 
an augmented "Nambu-Gorkov" spinor 




(15) 



where is the charge conjugated spinor. This doubling 
of the degrees of freedom requires careful attention to 
avoid double counting. 

To simplify the presentation of our method in this sec- 
tion, we shall ignore this complication and assume that 
^ contains a single set of operators with no duplicate 
degrees of freedom. We also consider only a single inter- 
action term, and subsume the momentum structure into 
the matrix structure. Explicitly dealing with these com- 
plications is straightforward and the details are presented 
in 57].^ 



^ To derive the full equations, one introduces the augmented struc- 
ture while imposing constraints on the matrices throughout the 
variation. For example, one must ensure that H = — CH"^C 
where = C^. To derive the proper momentum structure, 
one simply attaches momentum indices: for homogeneous states 
all correlations have the form (^t^^) <x <5'^(p — q) and the mo- 
mentum structure follows trivially. 
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We now express the variational Hamiltonian as 

iJo = ^^IH + S)* = *tHo*. (16) 

The matrix of variational parameters S may be thought 
of as the self-energy corrections. All of the two-point 
correlations are determined from the corrected Hamilto- 
nian matrix Hg, with the "anomalous" correlations being 
found off the diagonal: 

At finite temperatures, there is a one-to-one relationship 
between E and the matrix F~ = 1 — F'^. Armed with 
this result, the variational bound takes the explicit form^ 

n<- Trln[F-] - Tr[SF+] + 

P 

+ g (Tr[r^F+] Tr[rF+] + Tr[rtF-rF+]) . (17) 

From this, one may find the stationary points by by dif- 
ferentiating with respect to F^. At finite temperature 
this is formally equivalent to finding the stationary points 
by varying with respect to S. Differentiating with re- 
spect to S is complicated by the functional dependence 
and the result is not expressible as a simple matrix equa- 
tion. The conditions dil/dF^ = yield the fully self- 
consistent Schwinger-Dyson equations which may be ex- 
pressed as: 

S = 5 (^pt Tr[rF+] + r Tr[r^F+] + F^F^F - FF+F^) . 

(18) 

In principle, one may derive analytic expressions for these 
equations, but, in practise, the matrices are 72 x 72 and 
computing F^ analytically — even when many approxi- 
mations are made — is quite tedious. Instead, we simply 
use these expressions numerically. The required diag- 
onalization is then efhciently performed using standard 
numerical linear algebra tools. The traces involved in- 
clude momentum integrals, but for homogeneous states 
these are one-dimensional and thus also quite efhcient. 
The biggest challenge is to solve simultaneously the equa- 
tions present in H18I) . This is done by first projecting out 
the the limited subspaces describe in Appendix ^ and 
then employing a multi-dimensional root-finder. 

Since the search space is large (~ 45 parameters for the 
CFLK° states), traditional root-finders are prohibitive 
because they recompute the Jacobian at each step. Here 
we use a modified Broyden algorithm [s^ [s^ to provide 



For a fully self-consistent analysis, we must include the term 
— g (^tvjit^ ^ij/v]/^. These correlations vanish in this simplified 
analysis, but are included whe n the full augmented structure I15i 
is considered as discussed in |57l . Note also that momentum 
integration is implicit in the matrix multiplication and traces. 



a secant-like update to the Jacobian requiring far fewer 
function evaluations per step of the algorithm. 

In many cases, the Schwinger-Dyson equation H18|l con- 
verges through simple iteration. With charge neutrality 
constraints, this is often no longer the case, but the Broy- 
den update is sufficient to restore convergence. 



B. CFL at TUs = 0. 

As an example, consider the parity even CFL state. 
The self-consistency conditions are fully closed when one 
includes four variational parameters. There are two gap 
parameters A3 and Ag corresponding to the diquark con- 
densate Q , one chemical potential correction S/j. b to the 
baryon chemical potential and an induced off-diagonal 
chemical potential fioct ■ The quadratic Hamiltonian Hll|) 
can thus be expressed 

where 

[Arfs, = A,e"^''eabk + A^iS^S^, + S^5^,), (19) 
S fj, ^ ^S^B + fJ'oct, (20) 

and 

[fJ'oct]"'"' f3b = Moct ^^[-^"^laaiA^l/Sfc " 8 [A°] [A°]^fc^ . 

Most of this structure is all well-known and discussed 
many times in the literature, however, there has been 
no mention of the parameter /Xoct because most analyses 
neglect the self-energy corrections. 

Neglecting the correction to the baryon chemical po- 
tential is reasonable since it has little physical signifi- 
cance: it simply enters as a Lagrange multiplier to estab- 
lish a finite density. As such, the effective common quark 
chemical potential 

N = |Mlf =l{^^B + S^ib) (21) 

is the relevant physical parameter defining the Fermi sur- 
face. To compare states in the grand ensemble, however, 
one must fix the bare rather than the effective chemical 
potentials. This is what we have done in our calculations. 
Numerically, we find that the corrections SfiB cause /i^ 
to vary by only a few percent as we vary the perturbation 
parameters fiy and rus- 

There is no bare parameter corresponding to ^oct- 
Thus, it is spontaneously induced and should be treated 
on the same footing as A. To see that such a parameter 
must exist, consider changing to the "octet" basis using 
the augmented Gell-Mann matrices 

V^A = 2[X%ai^ic.a) (22) 
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where A° = l/-\/6- In this basis, the ofF-diagonal con- 
densate becomes diagonal with one singlet parameter 
4A6 + 2 A3 and eight octet parameters Ag — A3: 

/4A6 + 2A3 \ 



\ A6-A3/ 

It is clear that in the CFL, the singlet channel decouples 
from the octet channel: there is no symmetry relating 
these and the two gap parameters are related by the nu- 
merical value of the coupling g. This decoupling is also 
present in the chemical potential corrections. One linear 
combination corresponds to the identity: this corrects the 
baryon chemical potential 5^b- The other is the induced 

Moct- 

Numerically, we calibrate our model with this CFL so- 
lution. In particular, we chose our parameters to repro- 
duce the results of i4fl| . We use a hard cutoff at A = 800 
MeV, and a coupling constant chosen so that, with an 
effective quark chemical potential of Hq — 500 MeV one 
has a physical gap in the spectrum of Aq = A3 — Ag = 25 
MeV. This fixes the following parameter values which we 
hold fixed for all of our calculations: 

A = 800 MeV, (24a) 
gK^ = 1.385, (24b) 
/is/S = 549.93 MeV. (24c) 

With these parameters fixed, the fully self-consistent 
mean-field CFL solution has the following variational pa- 
rameters: 

A3 = 25.6571 MeV, 
Ae = 0.6571 MeV, 
S^iB/i^ -49.93 MeV, 

= -0.03133 MeV. 

As first noted in 0] , and discussed in 60] , the parameter 
Ag is required to close the gap equation, but is small be- 
cause the sextet channel is repulsive. In weakly-coupled 
QCD, Ag is suppressed by an extra factor of the cou- 
pling. This effect is numerically captured in the NJL 
model. The parameter /Zoct is also required to close the 
gap equation when the Hartree-Fock terms are included. 
It is also numerically suppressed. Recent calculations of- 
ten omit Ag and ^Loct- we see that this is numerically 
justified. 

The physical gap in the spectrum also defines the crit- 
ical hypercharge chemical potential for the CFL/gCFL 
transition 

nl^^Ao^ 25.00 MeV. (26) 



C. CFL at /iY,ms / 

Once one introduces a strange quark mass, one must 
introduce additional parameters. A simple way to deter- 
mine which parameters are required is to add the mass, 
then compute the gap equation and see which entries in 
the self-energy matrix are non-zero. By doing this for a 
variety of random values of the parameters, one can de- 
termine the dimension of the subspace required to close 
the gap equation and introduce the required parameters. 

In the case of the CFL state with non-zero hypercharge 
chemical potential, one only needs to introduce the pa- 
rameters and fis to ensure gauge neutrality: As dis- 
cussed in Sec, nil A1 none of the other parameters change. 
To go beyond the transition into the gCFL phase, how- 
ever, or to extend the results to non-zero temperature, 
one must introduce additional parameters. These in- 
clude the perturbation , the gauge chemical potentials 
/i3, /ig, and /Xe required to enforce neutrality, as well as 
nine gap parameters A12, A23, A13, A45, Agy, Asg, An, 
A22 , and A33 that fully parametrize the triplet and sex- 
tet diquark condensates. (These latter nine parameters 
correspond to the parameters (pi, ipi and Ui defined in 
reference \39\.) The additional parameters are chemical 
potentials similar to /ioct which are induced by the gap 
equations. The full set of parameters in discussed in Ap- 
pendix^ 

Adding a strange quark mass is more complicated. 
First of all, we need to introduce additional Lorentz 
structure. For homogeneous and isotropic systems, there 
are eight possible relativistic structures: 

A = 1 (g) (5^ + 75 (K) Sfi^ - 7o ® Sm - foJ^ ® Sm^, 

B = 7c75 «) A -I- 7c «) As -I- 7o7c75 + 7o7c ® '^s- 

Introducing quark masses requires one to introduce the 
additional Lorentz structure k |60j | to close the gap equa- 
tions, but these are found to be small. In total, one re- 
quires about 20 parameters to fully parametrize the CFL 
in the presence of a strange quark mass (see Table Hvjl . 

With the inclusion of a bare quark mass TOj one in- 
duces a chiral condensate (V'V') which in turn generates a 
correction to the quark mass. The resulting parameter in 
Ho ((TT|l is the constituent quark mass Mg which appears 
in the dispersion relationships for the quarks. It is this 
constituent quark mass that must be used when calculat- 
ing the effective chemical potential shift l@J. Generally 
the constituent quark mass is quite a bit larger than the 
bare quark mass parameter m^. For example, close to 
the phase transition, we have « 83 MeV while the 
constituent quark mass is Mg ~ 150 MeV (see Table ITvll . 
We have checked that our calculations are quantitatively 
consistent with the calculations presented in |6lj in this 
regard. 
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D. CFLK° 

Applying a kaon rotation to the CFL state breaks the 
parity of the state, and mixes the parity even parameters 
m, A and k with their parity odd counterparts /X5, 
m5, A5 and K5. The full set of parameters and typical 
numerical values is presented in Appendix fXI 

V. LOW-ENERGY EFFECTIVE THEORY 

To describe the low-energy physics of these models, 
we follow a well established procedure: identify the low- 
energy degrees of freedom and their transformation prop- 
erties, identify the expansion parameters (power counting 
scheme), write down the most general action consistent 
with the symmetries and power counting, and determine 
the arbitrary coefficients by matching to experiment or 
another theory. In our case, we will match onto the mean- 
field approximation of the NJL model. The resulting low- 
energy effective theory has been well studied 18, 62]: we 
use this presentation to establish our conventions, and 
to contrast the effective theory of QCD with that of the 
microscopic NJL model. 

A. Degrees of Freedom 

The coset space in the NJL model is isomorphic to 
U(3)(g)U(3). This can be fully parametrized with two 
SU(3) matrices X and Y and two physical phases A and 
V which one can physically identify with the condensates: 



VytI[X],^ cx eabce^p^ hl'T^^'l) , (27a) 

yyUt[Y],^cxeafccea/37(V'?rV'^^>- (27b) 
These thus transform as follows: 

X^FlXC^, (28a) 

Y ^ FflYC^ (28b) 
A ^ eMOR~er.)ji^^ (28c) 

V ^ e^'^''''+''^^V. (28d) 
Note that the condensation pattern X = Y = 1, A = 



y = 1 is unbroken by the residual symmetry where = 
Fj^ = C and also by the Z2 symmetries where 0l,9ji — 
±7r. This is the reason for the extra factor of two in 
the phases. In QCD the degrees of freedom are similar, 
but one must consider only colour singlet objects. Thus, 
the low-energy theory for QCD should include only the 
colour singlet combination 

S = XY^ ^ FlSF^ (29) 

and the colour singlet phases A and V. Note also that 
these have the following transformation properties under 



parity 

X«->Y, A<->A^ S<->S^ (30a) 
The field content of the effective theories is thus: 

H, 7y': Two singlet fields corresponding to the U(l) 
phases of A and V. The field associated with y is a 
scalar boson associated with the superfluid baryon 
number condensation. We shall denote this field H. 

The field associated with ^4 is a pseudo-scalar boson 
associated with the axial baryon number symmetry 
and shall be identified with the rj' particle. As dis- 
cussed in Section^ the axial symmetry symmetry 
is anomalously broken in QCD and the 77' is not 
strictly massless due to instanton effects, but these 
are suppressed at high density. We ignore these ef- 
fects. Our NJL model thus contains no instanton 
vertex and our low-energy theory will contain no 
Wess-Zumino-Witten terms |63, It would be 
interesting to include both of these terms and re- 
peat this calculation as these effects are likely not 
small m. 

7r°: Eight pseudo-scalar mesons 7r° corresponding to the 
broken axial flavour generators. As colour singlets 
these remain as propagating degrees of freedom in 
both QCD and NJL models. These have the quan- 
tum numbers of pions, kaons and the eta and trans- 
form as an octet under the unbroken symmetry. 

0": Eight scalar bosons cj)'^ corresponding to the bro- 
ken coloured generators. These are eaten by the 
gauge bosons in QCD and are removed from the 
low-energy theory. This gives masses to eight of 
the gauge bosons and decouples them from low- 
energy physics. In the NJL model these bosons 
still remain as low-energy degrees of freedom, but 
decouple from the colour singlet physics when one 
properly enforces colour neutrality. 

There are additional fields and effects that should be con- 
sidered as part of a complete low-energy theory, but that 
we neglect: 

1. The appropriately "rotated electromagnetic field" 
associated with the unbroken U(l)g symmetry re- 
mains massless. Both the CFL and CFLK" states 
remain neutral with respect to this field, however, 
and we do not explicitly include it in our formula- 
tion. 

2. The leptons are not strictly massless, but the elec- 
tron and muon are light enough to consider in the 
low-energy physics. In particular, they contribute 
to the charge density in the presence of an electric 
chemical potential and at finite temperature. In 
this paper, leptonic excitations play no role since 
we consider only T = and both CFL and CFLK° 
quark matter is electrically neutral for /ie = 0. The 
leptons play an implicit role in fixing /ig such that 

= in both insulating phases. 
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To be explicit, we relate all of the dimensional physical 
fields H, 7]', (ff and tt" to the phase angles through their 
decay constants; H = JrH, t]' = fr^'V' , 0° = Z^^", and 
7r° — /ttTt". The two U(l) phases angles have a slightly 
different normalization because of the normalization of 
the generators. This normalization is chosen to match 
the kinetic terms in the original theory and matches 
113: fj' = VeiOR - 61l), and H' = V6{eR + 9l). The 
realization of these transformations in the microscopic 
theory is 



Kinetic Terms 



ip exp ■ 



-HI - 77^75 
2^6 



vr 



V' (31) 



where 



fa 



(i± 75)® (-a:)® 1/2, 

1 (g) 1 (g) Aa, 

f?i - f2 = 75 ® i-K) ® 1, 

(f^-c'^)/2 



(32a) 
(32b) 
(32c) 
(32d) 
(32e) 



X 


= exp{ 


— ztt^Aq} exp 


Y 


— exp { 


171-" Aa} exp 1 


A 


= cxp| 




V 


= exp| 




S 


= exp{ 





and the corresponding realization in the effective theory 

is 



(33a) 

(33b) 

(33c) 

(33d) 
(33e) 



B. Power Counting 

In addition to Aqcd which separates the three light 
quarks from the heavy quarks, there are two primary 
scales in high-density QCD: the quark chemical poten- 
tial fiq and the gap A. In the NJL model there is also a 
cutoff and the coupling constant: these are related by the 
gap equation when one holds /z and A fixed and the qual- 
itative physics is not extremely sensitive to the remaining 
renormalization parameter. 

Our low-energy theory is an expansion in the en- 
ergy/momentum of the Goldstone fields. Thus, the ex- 
pansion is in powers of the derivatives with respect to 
the scales /i and A. In this paper, we shall only consider 
leading order terms: Systematic expansions have been 
discussed elsewhere (see for example (fi^). 



To construct the low-energy theory we follow ^3 ^-^d 
introduce coloured currents 



X (v) 



''x ^ ''y 



In the presence of a finite density, we no longer have 
manifest Lorentz invariance and must allow for additional 
constants into our spatial derivatives 



to account for the differing speeds of sound. This pa- 
per will be concerned with static properties, so we can 
neglect these. In principle, one must also match these co- 
efficients V. In QCD this matching, along with other coef- 
ficients, has been made with perturbative calculations at 
asymptotic densities [23, H^l . Our theory and states still 
maintain rotational invariance. Thus, to lowest-order we 
have m 



4 

, 3/i 



A- 



)^^%„)^^-^Tr[JtJ 



J TT 

4 

4 



Tr[J^lJ_^]- 



H 



+ 



The neglected terms are of higher order in the derivative 
expansion. Note that our normalizations have been cho- 
sen so that this expression is canonically normalized to 
quadratic order in terms of the dimensionful fields. 

The division of £eff is natural 0| because it separates 
out the colour singlets. £s depends only on S for exam- 
ple: 



Thus, with the exceptions noted above, the lowest-order 
low-energy effective theory of massless Nf = 3 QCD is 



C 



QCD 



H 



r 



(34) 



whereas the NJL model proper must also include £4 



D. Perturbations 

We shall now consider two types of perturbations: 
chemical potentials and quark masses. To deal with these 
perturbations, we note that they enter the microscopic 
Lagrangian as 

£sB = V'Im^^l + i'Rfj'^i'R + i'RMiPL + V-ImVa- 
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These terms break the original symmetries of the the- 
ory, but one can restore these symmetries by imparting 
the following spurion transformations to the masses and 
chemical potentials 



M 

..L 



R 



(FI< 



. C)M(F* d 
C)/x^(F2 c 
'C)/x«(F^,, 



)C)t, 
5C)t. 



(35a) 
(35b) 
(35c) 



The transformation M — M preserves the residual 
Z2 symmetries. This prevents odd powers of the mass 
terms from appearing in the chiral effective theory. In 
particular, the linear term dominant in the vacuum is 
forbidden, resulting in an inverse mass-ordering of the 
mesons [23, 123 with the kaon being the lightest particle 
at high density. 

All these symmetries must be restored in the effective 
theory: we are only allowed to couple these parameters 
to the fields in ways that preserve the global symmetries. 
To lowest-order, this greatly limits the possible terms in 
the effective theory. 

In the case of the chemical potentials, we can go one 
step further by noting that the perturbations always ap- 
pears in combination with the time derivative 



(36) 



One can thus promote the chemical potentials to a tem- 
poral component of a spurion gauge field and render the 
symmetries local in time: 



/X ^ (F* ® C) (/X + ido) (F* ® C)t. 



(37) 



The effective theory must also maintain these local sym- 
metries. One concludes that the chemical potential per- 
turbations can only appear through the introduction of 
covariant derivatives in the effective theory. In particular, 
consider adding independent colour and flavour chemical 
potential terms: 



L.R 



(38) 



where jip and /Xp are traceless 3x3 matrices. From these 
we may construct the following quantities that transform 
covariantly: 



VoX = 


9oX 




*X4 




(39a) 


VoY = 


9oY 


■f 


*Y4 




(39b) 


VoS = 










(39c) 




(9o + 




V, 




(39d) 


VoA = 


(9o + 


2ifiA) 


A, 




(39e) 



where /iy = + is a small adjustment of the baryon 
chemical potential ^_b/3 and = M-R ~ Ml is the "axial 
baryon" chemical potential. For the rest of this paper, 
we shall only consider vector chemical potentials that are 
real and symmetric: /x^'^ — — ■■* — 



these restrictions, the static potential in the effective the- 
ory is 



■Tr 



(40) 



to lowest-order. The terms omitted include terms of 
higher order in the perturbation and small corrections 
due to the explicit violation of the "local" spurion sym- 
metries by the cutoff. 



E. Charge Neutrality 

As discussed in 0, the gauge-invariance of 
QCD implies that homogeneous states must be colour 
neutral. Non-zero static colour sources A^(p = 0) cancel 
the tadpole diagrams ensuring neutrality. These sources 
enter the NJL calculation as Lagrange multipliers to en- 
force neutrality. 

One can see explicitly how these arise in the context of 
the effective theory. The gauge fields effect the local sym- 
metry and thus couple through the derivatives in exactly 
the same way as the spurion coloured chemical poten- 
tials: oc Enforcing gauge-invariance induces 
an effective coloured chemical potential that makes (|40f) 
stationary with respect to variations of the gauge field, 
and thus equivalcntly,with respect to traceless variations 
of /ip. Thus, we sec that, to lowest-order [Tsl Is^. Issf 



= -i (xVfX + yVfY 



(41) 



Inserting this into the ()40() , and considering only traceless 
perturbations, we see that the colour dependence drops 
out of the effective theory and we are left with the static 
effective potential involving only the colour singlet fields: 



(42) 



Mf = Mf = Mf stc. With 



In order to reproduce the physics of this in the NJL 
model, however, we must remove the coloured degrees 
of freedom. This is done by introducing colour chemical 
potentials to the NJL model as Lagrange multipliers and 
using them to impose colour neutrality 45, 54ii.l55l]. This 
removes the colour dependence in the NJL model to all 
orders in the same way as it removes the colour depen- 
dence in (|41|l to lowest-order. (In general, it is not suffi- 
cient to impose colour neutrality: one must also project 
onto colour singlet states (as well as states of definite 
baryon number). This projection is important for small 
systems, but likely has negligible cost for thermodynam- 
ically large systems such as neutron stars. See for an 
explicit demonstration of this in the two-flavour case.) 

The quarks also couple to the photon, and so we also 
must enforce electric neutrality. Enforcing electromag- 
netic gauge-invariance will likewise induce an electric 
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FIG. 7: Numerical difference in energy densities between 
the kaon-condensed CFLK*^ state and the CFL state at fi- 
nite hypercharge potential /iy obtained from our microscopic 
NJL calculation. The units are scaled in terms of the quark 
chemical potential /i, = 500 MeV and the critical hypercharge 
chemical potential ^.y = 25 MeV. The quantities plotted were 
chosen so that the relationship will be linear if our calcula- 
tion agrees with the effective theory result l|45|l . The slope of 



the line is m 



-0.018 from which we can deter- 



mine the effective theory parameter /,r ~ O.lO/i,. This is in 
good numerical agreement with the perturbative QCD result 



U ~ 0.209At5 



The dashed extension shows the com- 



parison between the CFLK° potential and the CFL potential, 
but beyond 1.0, the CFL becomes the gCFL and the energy 
dependence changes. We have not calculated the gCFL po- 
tential in this paper, but plot this extension to emphasis that 
the CFLK'' persists beyond the CFL/gCFL transition point 
at 1.0. 



chemical potential /.te that ensures electric neutrality. It 
turns out that both the CFL and the CFLK° quark mat- 
ter are neutral under a residual charge Q (both are Q 
insulators). This means that one has some freedom in 
choosing the chemical potentials used to enforce neutral- 
ity. In particular, prior to the onset of gapless modes, 
one may choose these combinations so that = 0. This 
is naturally enforced by including charged leptons in the 
calculation. 

Once a charged excitation becomes gapless, the mate- 
rial becomes a conductor and a non-zero /^e is required 
to enforce neutrality. The phase transition to the gCFL 
and gCFLK" is defined by exactly such a charged excita- 
tion. In this paper, we shall only consider the insulating 
phases, and thus simply set /Xg = 0. For further discus- 
sions of the metal/insulator properties of the CFL and 
gCFL we refer the reader to |37Ll39ll40|. 



F. Kaon Condensation 

We are now in a position to argue for the existence of a 
kaon-condensed state. Consider performing an axial K*^ 
rotation on the parity-even CFL state. This is effected 
using (|31|l in the microscopic theory and using H33|l in 
the effective theory with the parameter tt^ = 9. Such a 



FIG. 8: Numerical difference in energy densities between 
the kaon-condensed CFLK" state and the CFL state at finite 
strange quark mass nis obtained from our microscopic NJL 
calculation. The units are scaled in terms of the renormal- 
ized quark chemical potential /i, ~ 500 MeV and the crit- 
ical hypercharge chemical potential /ly = 25 MeV to facil- 
itate comparison with Figure |7| and to emphasize the lin- 
ear relationship implicit in 147^ . The slope of the line is 
m ~ — /^/2/iq ~ —0.028 which gives an effective « 0.21/iq 
which is consistent with our previous results. In comparison 
with Figure |7| the CFL— >gCFL transition occurs somewhat 
earlier because the gap parameters are reduced with increas- 
ing strange quark mass. The curve cannot be extended as in 
Figure |7| because the free-energy of the CFL is no longer a 
constant as it was with a hypercharge perturbation. 



state is now described by 



/l 

cos(0) isin(0) 
ism[e) cos{e) 



(43) 



In the presence of a hyperc harge c hemical potential, the 
effective potential becomes IMl 



Vi9) 



W(^)-l) 



(44) 



We see that this has a minimum for = ±7r/2: this is 
the state with maximal K° condensation. We can also di- 
rectly compute the difference in the thermodynamic po- 
tential densities between the CFL state and the CFLK° 
state: 



CFLK" 



CFL 



(45) 



Armed with this relationship, we can now turn to the 
microscopic calculation and determine the coefficient 
In Figure we plot our numerical results so that the 
linear relationship (|45|l is evident. From the slope of the 
relationship we find that 



(46) 



We note that this is in good numerical agreement with 
the perturbative QCD result HI IS of = 0.209^^,. 
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FIG. 9: Chemical potentials required by the NJL model to 
enforce colour neutrality in the CFLK" phase with finite hy- 
percharge chemical potentials. Note that the effective theory 
relationship 1501 is satisfied from small chemical potentials. 
The linear deviation seen here reflects the missing terms in 
the effective theory that are of higher order in the perturba- 
tion /iy with a linear deviation here corresponds to /.ty terms 
missing in l|50|l . 

This striking agreement between two very different mod- 
els arises from the fact that this coefficient is not very 
sensitive to the effects of the cutoff (which is different in 
the two theories) and gives encouraging support to the 
use of the NJL model to study QCD. 

The equivalent relationship in the case of a strange 
quark mass requires one to include mass terms in the 
effective theory (see for example 113), but the leading 
order effect can be determined by using the "effective" 
strange quark chemical potential iiy ~ ^s/C^l^q) that 
follows from ((SJ: 

W-^^CPL = -^(^) . (47) 

It is important here to note, however, that the strange 
quark mass affects the solution in such a way that the gap 
parameters change and self-energy corrections modify the 
quark chemical potential and the constituent quark mass. 
It is the renormalized parameters that appear in this re- 
lation and in the perturbative QCD result. Thus, as a 
function of the bare parameter Ws, we have Ms oc to^ 
and ^iq — fx m^. Thus, we should see a linear re- 
lationship between r^cFLK" ~ ^cfl and M^- We plot 
this relationship in Figure |H1 and extract the slope which 
gives the relationship w 0.21^g. This is in qualitative 
agreement with our previous result. The slight numerical 
disagreement is due to effects of the strange quark mass 
that are not captured by the chemical potential shift JJJ . 

We pause here to point out a discrepancy between 
our results and similar work by Buballa |67|. Our re- 
sults shown in Figure |H1 suggests that kaon condensation 
occurs for all values of in this simple model with 
rUu = TUd = whereas Buballa finds that kaon conden- 
sation is only favoured for to^ sufficiently large. If the 
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FIG. 10: Hypercharge density of the CFLK state in the 
presence of a hypercharge chemical potential /iy as obtained 
from our microscopic NJL model calculation. The units are 
scaled as in Figure |7| so that the relationship will be linear 
if the NJL model calculation agrees with the effective theory 
prediction 1511 1. By determining the slope of this relationship 
we have another way of determining the coefficient in the 
effective theory. The slope is f^/fiq « 0.037 which agrees 
with our previous determination of /,r ~ 0.19/ig. 



chemical potential shift were the only effect of a strange 
quark mass, then this would be inconsistent with H40|l. 
This is not, however, the correct expansion. Instead, one 
has, for maximal kaon condensation and m„ — rud = 0: 

^^CFLK" - f^CFL ^ -(4a6 + ag^ml - cm^ + O(m^). (48) 

A proper discussion of this is beyond the scope of this pa- 
per, but will be discussed elsewhere The term with 
coefficient c should be identified with the leading order 
contribution from ^ ; the term with coefficient ag arises 
from the sextet gap contribution Ag and is discussed 
in ^3 ; the term with coefficient ag is of higher order in 
perturbative QCD and so is usually neglected p^l23ll30l |. 

Thus, Buballa's results are consistent with the effective 
theory. The discrepancy is due to a different choice of 
parameters A ~ 100 MeV and A = 600 MeV compared 
with our parameters A 25 MeV and A = 800 MeV. 
With the large gap, one is further from the perturbative 
QCD regime and the quadratic term appears to play a 
significant role. In our analysis, the term — (4a6 + a^)m1 
is small compared with the cm^ term. Using Buballa's 
parameters, however, we qualitatively reproduce his re- 
sults. A further discussion of these effects will be pre- 
sented shortly ^68j. 

There are a couple of other consequences that follow 
directly from the effective theory. One is the value of the 
coloured chemical potentials required to enforce neutral- 
ity. In our microscopic model, we have fixed the gauge 
(unitary gauge) by setting X = = \/S for the axial 
rotations. The CFL state has X = Y = 1 while the 
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CFLK° state has 




(49) 



From H41I) we have the following relationships required to 
enforce neutrality 



= 0, 



(CFL), 



(50a) 



^,s = -l^,Y, M3 = -^Mi-, (CFLKO). (50b) 

We plot these relationships in Figure |21 Note that they 
only hold for small perturbations where the effective the- 
ory is valid: this plot also demonstrates a departure from 
the lowest-order effective theory as the perturbation is 
increased. 

As a final demonstration of the effective theory, we 
calculate the hypercharge density. This is obtained by 
varying the thermodynamic potential with respect to the 
hypercharge chemical potential: 



ny 



~Jl^ly{cos\e)~l). 



(51) 



There should be no hypercharge density in the CFL state 
and a density of ny = f^fxy in the CFLK" state. In- 
deed, the CFL supports no hypercharge density with 
nu = rid = ris- The hypercharge density of the CFLK'^ 
phase is shown in Figure lTUI and provides another method 
of extracting — 0.19^,. 



G. A Note on the Meaning of 1/(61) 

We make a few remarks here about the meaning of 
the effective potential V{9). In particular, one might be 
tempted to try and compute the functional form of V{9) 
in the microscopic theory to facilitate matching with the 
effective theory. Such an approach will generally fail be- 
cause one is allowed to pick an arbitrary parametrization 
of the Goldstone fields as long as they leave the kinetic 
terms unaltered [s^ [t^ . Physical quantities must be in- 
variant under this change of parametrization: thus the 
spectrum about the minimum, densities, and energy dif- 
ferences are reasonable quantities to compare in each the- 
ory. The general form of the effective potential away from 
the stationary points, however, is rather arbitrary. 

As an example: consider starting with the parity even 
CFL state in the presence of a finite fiy. This state 
corresponds to a stationary point of the effective poten- 
tial and is a self-consistent solution to the gap equations. 
One can then form a continuum of "kaon rotated" states 
\9) by applying the broken symmetry generators to this 
state. One might expect to find V{6) by computing the 
energy of these states, but instead one finds an expres- 
sion that is only valid locally about the stationary point. 
The reason is twofold: First, there is not a unique "kaon 



rotated" state \9). This state has many other parame- 
ters corresponding to other "directions" (such as the gap 
parameters A, the chemical potential corrections etc.) 
The only way to uniquely determine these is to solve 
the gap equations, and these only have well-defined so- 
lutions at stationary points. Second, the generators of 
the pseudo-Goldstone bosons in the presence of pertur- 
bations are not the same as the generators of the true 
Goldstone bosons in the unbroken model: the pseudo- 
Goldstone bosons have some admixture of these other 
"directions" . 

This becomes even more evident when you analyze the 
CFLK'' state with a large perturbation: one can try to 
"undo" the kaon rotation by applying the appropriate 
symmetry generators to minimize the parity violating 
condensates, but one finds that there is no way to do 
this. One must also transform the other parameters in 
order convert a CFLK" state back to a parity even CFL 
state. 



VI. CONCLUSION 

We have explicitly found self-consistent solutions 
within a microscopic NJL model exhibiting the feature of 
kaon condensation in a colour-flavour-locked state. Us- 
ing these solutions, we have demonstrated that by prop- 
erly enforcing gauge neutrality, one can remove the extra- 
neous coloured degrees of freedom from the NJL model 
and effectively model kaon condensation in high-density 
QCD. In particular, the microscopic calculations can be 
matched onto the low-energy effective theory of QCD. 
We determined /jr = 0.19^^ which is in good numerical 
agreement with the perturbative QCD result. 

Furthermore, our solutions are fully self-consistent: no 
approximations have been made beyond the mean-fleld 
approximation and restricting our attention to isotropic 
and homogeneous states. We find that our results 
agree qualitatively with both the expected properties 
of the CFLK° phase based on effective theory calcula- 
tions, and with the previous numerical calculations of 
the CFL/gCFL transition. 

Quantitatively we find that the phase transitions oc- 
cur at slightly smaller parameter values than previously 
found in the literature. Concerning the CFL/gCFL tran- 
sition, we find that the gap parameters are reduced by a 
few percent compared with those presented in and 
subsequently, the critical Mg is also a few percent lower. 
Concerning the CFLK"/gCFLK° transition, we find that 
the transition occurs about a factor of 1.2 higher than the 
CFL/gCFL transition. This is in qualitative agreement 
but quantitative disagreement with the factor of 4/3 cal- 
culated in ^34] . 

The next step is to use this microscopic model to de- 
termine the phase structure of high-density QCD in the 
region where the gapless modes appear. We suspect that 
the gCFLK° state will survive somewhat longer than the 
gCFL state on account of its lower condensation energy. 
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but a quantitative comparison is required. Extrapola- 
tion to finite temperature is also a trivial extension in 
our formalism. 

A somewhat more challenging direction is to consider 
the effects of instantons and finite up and down quark 
masses and investigate other forms of meson condensa- 
tion. Preliminary investigations indicate, however, that 
the number of parameters required to close the gap equa- 
tions in the presence of arbitrary meson rotations may be 
prohibitively large to continue with fully self-consistent 
calculations. This should still be tractable with carefully 
made approximations. 
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APPENDIX A: FULL PARAMETRIZATION 

In this appendix, we give the full parametrization used 
to analyze the K° condensed states. First, we must in- 
troduce a full set of diagonal chemical potentials. One 
approach would be to introduce the 9 individual quark 
chemical potentials, but certain linear combinations cou- 
ple to relevant physics. We fix the overall density by 
fixing the baryon chemical potential fis- Then we must 
enforce gauge neutrality, so we introduce /ig which cou- 
ples to the electromagnetic field, and the diagonal colour 
chemical potentials fi^ and /ig. The rest of the chemi- 
cal potentials are chosen to be orthogonal to these. Here 
then are the diagonal elements of the diagonal chemical 
potentials expressed as tensor products of the flavour and 
colour structure : 
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An alternative set of chemical potentials includes the hy- 
percharge chemical potential fiy instead oi fif. These are 
no longer orthogonal, but are still linearly independent. 
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The diagonal mass corrections (chiral condensates) do 
not couple to any external physics, so we simply use the 
nine quark mass corrections ((5m„r, ^rriug, Srriub, Srridr, 
dnidg, Srudb, Sm^r, Snisg, Snisb)- 

The rest of the parameters are described in the fol- 
lowing matrices. These appear more condensed when ex- 
pressed in the basis described in [i^l where the quarks are 
ordered (ru, gd, bs, rd, gu, rs, bu, gs, bd). In this basis, 
the matrices corresponding to the variational parameters 



A = 1 ® (5/x 75 ® (5/^5 - 7o ® (5m - 7575 ® (Jms, 

B = 7c75 «) A -f 7c. «) A5 -I- 7o7c75 ®'^ + 7o7c ® «^5- 
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In order to allow for a computer to enumerate the pa- 
rameters, we introduce a systematic method for labelling 
the parameters. First, we use one of the names /x, /i^, 
m, 771^, A, A^, K, or corresponding to the structure 
given above. We then use a two-digit index to specify 
which elements are non-zero and an i indicates that the 
specified element is i rather than simply 1 . The symmet- 
ric entry must also be set so that the resulting matrix 
'y ® II is either Hermitian or anti-symmetric depending 
on whether or not it parametrizes A or B respectively. In 
total, there are 666 independent matrices. For example 

/O 1 --A 
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The reason that and behave differently than the 
others is that, while 1, 75, and 7o Hermitian, 7075 
anti-Hermitian. Likewise, while 'Jc'YbJ ^'^'^ ToTcTs 
are anti-symmetric, 7o7c is symmetric. Again, recall 
that these are all specified in the "unlocking" basis which 
is ordered as 



ru, gd, bs, rd, gu, rs, bu, gs, bd. 



(A5) 



The parity even CFL state with no mass or hypercharge 
is expressed in terms of this parametrization as 
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Here are some comparisons with other conventions in the 
literature. Alford, Kouvaris, and Rajagopal j^J] intro- 
duce Ai, A2 and A3_ which are all related to the attrac- 
tive anti-symmetric 3 channel: 



(A7a) 
(A7b) 
(A7c) 



Riister, Shovkovy, and Rischke introduces the pa- 
rameters 0, if and a which include the repulsive sym- 
metric 6 channel parameters: 
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Finally, Buballa [63 uses only the following parameters 
to parametrize the meson condensed phases: 
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In Tables iHl lmlUvl andE| we give the numerical values 
of the parameters for each of the states displayed in Fig- 
ures 13 and respectively. We only list the non-zero 
parameters: the other parameters are zero. 
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+25.657 


+ 


+25.657 


As 





+0.65709 





+0.65709 



TABLE IL Parameters required for a self-consistent parity- 
even CFL solution in the presence of a hypercharge chemical 
potential. These values correspond to the dispersions shown 
in Figure Q All values are in MeV. The first column labelled 
"Bare" gives the fixed bare parameters that enter the Hamil- 
tonian Q. The column labelled "Correction" is the contri- 
bution from the self-energy. The sum of the columns is the 
value that enters the quadratic Hamiltonian IIH . 
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My = 


-- 0.50My 


Mi- 


= 1.20My 




Par am. 


Bare 


Correction 


Bare 


Correction 






/is/3 


+549.93 


-49.932 


+549.93 


-49.947 






Hy 


+12.5 


-1.0687 


+30 


-2.5931 






Me 





+0.53436 





+1.2965 






M3 


-6.4772 


-0.00000 


-16.346 


-3.331 X 10^^ 






MS 


-3.2386 


-0.00000 


-8.1729 


-1.665 X 10"^ 






Me3 





+0.02421 





+0.027856 






Me8 





+0.0080699 





+0.0092853 






M/3 





+0.035998 





+0.084801 






M/8 





+0.011999 





+0.028267 




/il2 


= Ml8i 





+0.016852 





-0.047357 






— Ml9i 





+0.11902 





+0.19694 




M23 


= M89 





+0.046967 





+0.046617 











+0.088616 





+0.14526 






— M38i 





+0.046967 





+0.046617 






M39i 





+0.0030234 





-0.065811 











+0.11479 





+0.27514 






All 





+0.64468 





+0.5851 


A22 = 


-Ass = 


-Ais, 





+0.32265 





+0.31566 


A33 = 


-A99 = 


— ^39i 





+0.33523 





+0.34523 




A12 = 


-Afs, 





+9.6383 





+10.138 




Al3 = 


-Afg, 





+8.9893 





+8.5746 




A23 = 


-Agg 





+12.789 





+12.605 




A45 = 


A^ 





-9.1811 





-9.7012 




A67 = 


A5 
— ^47i 





-8.5229 





-8.1128 




^29i - 


- A^ 

- ^38i 





-0.33242 





-0.35025 



TABLE III: Parameters required for a self-consistent CFLKO solution in the presence of a hypercharge chemical potential. 
These values correspond to the dispersions shown in Figure |21 All values are in MeV. The first column labelled "Bare" gives 
the fixed bare parameters that enter the Hamiltonian (|7J. The column labelled "Correction" is the contribution from the 
self-energy. The sum of the columns is the value that enters the quadratic Hamiltonian llllll . 
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M^/{2fi) = 0.83^y 


Param. 


Bare 


Correction 


Bare 


Correction 


/is/3 


+549.93 


-48.952 


+549.93 


-48.316 


A^s 


-12.649 


-0.00000 


-20.95 


+1.469 X 10"^ 


Me3 





+0.022617 





+0.022026 


Me8 





+0.0074542 





+0.0072087 


M/3 





-0.022617 





-0.022026 


M/8 





+0.022362 





+0.021626 


Ml2 





+0.090467 





+0.088103 


Ml3 = M23 





+0.088963 





+0.085669 







+0.15778 





+0.19366 


mug = 





+0.17255 





+0.2117 







+0.15604 





+0.19155 




+61.843 


+50.029 


+80 


+64.267 


rusb 


+61.843 


+50.079 


+80 


+64.329 







-0.014765 





-0.018037 


mi3 = 77123 





+0.026496 





+0.032723 


All = A22 





-0.62077 





-0.59745 


A33 





-0.64043 





-0.62926 


A12 





-12.914 





-12.753 


A13 — A23 





-12.639 





-12.302 


A45 





+12.293 





+12.155 


AgT = Agg 





+12.011 





+11.693 


TCll 





+3.8762 X lO"** 





+5.1011 X 10"'^ 


K22 





+3.8762 X 10^** 





+5.0102 X 10"*^ 


K33 





+0.078913 





+0.098773 







+0.0017234 





+0.0020923 


Kl3 — K23 





+0.52751 





+0.66061 


«:45 





-0.0017195 





-0.0020872 


K67 = K89 





-0.48835 





-0.61184 



TABLE IV: Parameters required for a self-consistent parity even CFL solution in the presence of a strange quark mass. These 
values correspond to the dispersions shown in Figure 13 All values are in MeV. The first column labelled "Bare" gives the fixed 
bare parameters that enter the Hamiltonian (|7|l. The column labelled "Correction" is the contribution from the self-energy. 
The sum of the columns is the value that enters the quadratic Hamiltonian IIH . For example, the right set of data (just slightly 
before the CFL/gCFL transition) has a bare (current) strange quark mass of 80 MeV. This corresponds to a constituent quark 
mass of 80 + 64 « 144 MeV. (Note that there is a slight difference for the blue constituent quark masses because of the presence 
of the coloured chemical potential ns required to enforce neutrality.) 
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M^/{2fi) = O.SO^y- 


M^/{2fi) = 0.84^y 


Param. 


Bare 


Correction 


Bare 


Correction 


/is/3 


+549.93 


-48.951 


+549.93 


-48.09 




-6.6057 


-0.00000 


-13.002 


-8.7 X 10"** 


MS 


-3.3029 


-0.00000 


-6.5008 


-4.35 X 10~* 







-0.53978 





-1.0238 







+0.023555 





+0.025582 







+0.0078516 





+0.0085274 


M/3 





+0.035852 





+0.066184 


M/8 





+0.011951 





+0.022061 







+0.014623 





-0.026924 


Ml3 = ^J■19^ 





+0.11615 





+0.16265 


M23 = M89 





+0.044694 





+0.042908 


/^28i 





+0.086431 





+0.11903 


5 5 
M29i — M38i 





+0.044644 





+0.042927 







-0.001146 





-0.047235 







+0.11528 





+0.21622 


f^ug — fTT'uh 





+0.0077476 





+0.0093125 


rUdr 





+0.17093 





+0.21503 


rudg — rudb 





+0.15529 





+0.19536 




+61.637 


+50.279 


+85 


+69.188 


rrisg — rrisb 


+61.637 


+50.335 


+85 


+69.301 


mi2 = -mi8i 





-0.010709 





-0.013162 


mi3 = -mfgi 





+0.029548 





+0.053931 


^123 = ?Tl89 





+0.0134 





+0.017695 


m|g, = -mlsi 





-0.012525 





-0.01659 







-0.0077707 





-0.009337 


All 





+0.60701 





+0.53896 


A22 = —Ass 





+0.30409 





+0.28263 


A33 = — Agg 





+0.32668 





+0.32494 


A12 = -Afsi 





+9.4718 





+9.6392 


Ai3 = -Afg, 





+8.614 





+8.0104 


A23 = — Asg 





+12.282 





+11.726 


A45 = — Algi 





-9.0411 





-9.2445 


AfiT = — A47i 





-8.1684 





-7.5847 


-^281 





-0.30382 





-0.28208 


^2gi — ^38i 





-0.31799 





-0.31604 


Aig, 





-0.32038 





-0.31321 


K22 = — K88 





-1.7029 X 10"^ 





-1.3167 X 10~^ 


K33 = — Kgg 





-0.040639 





-0.055178 


K12 = Acfgi 





-0.00066095 





-0.00090186 


«;i3 = Kl9i 





-0.36095 





-0.46324 


K23 = —^89 





-0.50626 





-0.66374 


5 

'1^45 — —l^hdi 





+0.00075037 





+0.001152 


K67 = K47i 





+0.33313 





+0.4268 


'tigi = — '*38i 





-0.019786 





-0.026852 



TABLE V: Parameters required for a self-consistent CFLKO solution in the presence of a strange quark mass. These values 
correspond to the dispersions shown in Figure |S| All values are in MeV. 



